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Abstract 

The naive use of higher order perturbation theory leads the left-hand cut 
O ' integrals in -kit dispersion relations divergent. This problem is discussed 

and solved. Also we point out that the Adler zero condition imposes three 
constraints on the dispersion relations. The a pole position is determined 



using the improved method, M a = 483 ± 13MeV,r CT = 705 ± 50MeV. The 
H-s . scattering length parameter is found to be in excellent agreement with the 

' experimental result. 

<n : 

^ : PACS numbers: 11.55.Bq, 14.40. Cs, 12.39.Fe 

>• \ Key words: nn scattering, dispersion relation, a meson 

' In Ref. PIE], we have established a new dispersion representation for the partial 

wave tttt scattering S matrix. The key point for setting up the dispersion represen- 
ts! . tation is the observation that both the real part and the imaginary part of S denned 
in the physical region are analytic functions on the cut plane expressed in terms of 
poles, the left hand cut integrals and the kinematic factor, p = yjl — Am^/s. To be 
specific, the real part of S, cos(25,r) (where 5 n is the scattering phase shift) and the 
imaginary part of the S matrix, sin(25 7r ) satisfy the following dispersion relations: 

> i sin(2<y = pF , 

><■ F(s) = a + T - V - 

S i i msi){s - 8i) f 2ip(z I j I )S'(z I /)(s - zf) 



1 f lm L F{s') 



7T JL S 1 



ds , + i_njn^m ds , 

7T JR. S' — S 



and 



cos^) = F = a + J2 



2{s- 8i ) 



, \- 1 , 1 flmJVl 1 f lm R F(s') , 

where Si denote the possible bound state pole positions and $ are the corresponding 
residues of S; zj 1 denote the possible resonance pole positions on the second sheet. 
The integrals denote the cut contributions, L=(— oo,0] is the left hand cut (l.h.c.) 
and R starts from the A A threshold once the Air cut are neglected, a and a denote 
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the subtraction constants, and one subtraction to the integrals in above expressions 
is understood. With these dispersion relations we can then generalize the single 
channel unitarity relation S + S = 1, which is only valid in the single channel physical 
region when s is real, to the whole complex s plane PQ: cos 2 (25 7r ) + sin 2 (25„-) = 1. 
The latter is equivalent to the well known generalized unitarity condition in quantum 
mechanics but was firstly discussed in field theory in Ref. PQ. 

The above method is valid for any partial wave scattering. It is however worth 
pointing out that in the scattering process with a non-vanishing angular momen- 
tum, J, restrictions among parameters should exist to ensure the threshold behavior, 
S n (s) oc k 2J+1 , where k = \/s — 4. In order to make use of Eqs. (JTJ) and in phe- 
nomenological discussions, a knowledge on the l.h.c. integrals is necessary. It is 
not very clear how to calculate these l.h.c. integrals in the nonperturbative scheme. 
Predictions on the left hand cuts from nonperturbative models like the Pade approx- 
imation are not always trustworthy Therefore results from chiral perturbation 
theory (CHPT) are used in estimating these integrals via the following formula, 

Im L F = 2Im L Re R T(s) = 2Im L T - 2Re L Im R T , (3) 

Im L F = -2p(s)Im L Im R T(s) , (4) 

since 

F = 2Re R T , F=l-2plm R T. (5) 

In Eqs. (J3J), (HJ) and ©, T is the partial wave scattering T matrix: S = 1 + 2ipT. 
The quantity Ihil-F is estimated at 0(p 4 ) in [2] and Ihil-F vanishes at this order. A 
question naturally arises as how close is the 0(p 4 ) results to the real situation. Since 
at 0(p 4 ) higher resonances do not contribute to Ihil-F and Ihil-F, it is argued in 
that the contribution from the t channel resonance exchange to the l.h.c. is very 
small at low energies. It is confirmed by a calculation to the tree level p exchange 
diagram [3], which indicates that the contribution to the l.h.c. integral is numerically 
very small up to s = lGeV 2 . Since resonance contributions saturate the 0(p A ) nn 
interaction chiral Lagrangian [3] at low energies and they contribute Ihil-F and 
Im L F at 0(p 6 ), the above discussions may suggest that the 0{p A ) results of Im L F 
and Im L -F are good approximations to the real situations at s < lGeV 2 despite 
of the ambiguity in choosing the cutoff parameter when estimating the integrals. 
However, even though high energy contributions to the left hand cut integrals may 
be small, the availability of the systematic use of perturbation theory in estimating 
Iull-F and ImL-F needs to be proved. This issue is not as trivial as it looks like at 
first glance. 

To have an understanding to the problem occurring in using higher order per- 
turbation theory results when estimating Ihil-F and Ihil-F, let us focus on Eq. (@J). 
Since Iull-F vanishes at 0(p A ), the leading order contribution is of 0(p 6 ), 

Im L Im R T(s) = 2pT 2 Im L Re / jT 4 

= 2pT 2 Im L T 4 - 2p 2 T| . (6) 

To obtain the second equation we have made use of the perturbative unitarity re- 
lation, IhirT 4 = pT 2 2 where (and hereafter) the subscripts denote the order of the 
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chiral expansions. Taking into account the results from CHPT we find that when s 
approaches 0_, Ihil-F behaves as 0((1/a/^s) 3 ) due to the presence of the kinematic 
factor. The integration in the left hand integral in Eq. (J2J) is therefore divergent 
when s' — > 0_. If higher order results are used the problem is getting worse since 
there will be higher powers of \/\f—s. The same situation occurs in Iull-F when 
using higher order results from ChPT. We will demonstrate in the following that this 
problem is only a deceptive artifact inherited from perturbation theory and can be 
corrected. In fact, from rather general considerations ImL-F(s) and lm.^F(s) should 
behave as 0(y/—s) and 0(l/y/—s), respectively, when s approaches zero. Hence the 
left hand integrals in Eqs. ((TJ) and (J2J) are well defined quantities. 

In order to understand the behavior of Ihil-F(s) when s — > 0, we should firstly 
understand the behavior of Tj(s) and Im^Tj(s) as s — > 0. Since s = is a branch 
point for Tj(s) we first let s approaches from the positive side along the real axis. 
From the partial wave projection formula 

1 f° It 
T I r(s) = ; — / dtPj(l + — r)r 7 (s,t,u) , u = 4ml-s-t, (7) 

we conclude that the partial wave amplitude is regular at s = 0+, since in the 
unphysical region s = 0, < t,u < 4m\ the full amplitude T 7 (s,t,w) contains 
no singularity. CHPT results are consistent with the conclusion obtained from the 
general analysis: actually in CHPT, at each order of the chiral expansion, the limit 
Tj/ 2n )(0 + ) exist as demonstrated by the explicit calculation for n = 1,2,3 and it is 
reasonable to assume it exists for arbitrary n. For the behavior of IhilTJ(s) near 
s = 0, we have 

1 + (— 1) I+J r4m£-s Of 

Im ^M = 32,(1- Ll) L ^P A l + 7 ^-;)l m T t \s,t);s<0, (8) 

from which we conclude IhilTJ(s) ~ (V~ s ) 3 as s ~^ 0- : one factor s comes from 
the integral interval in the above equation, another factor of y/—s comes from the 
threshold behavior of ImT/(s,t) near t = 4m\. This observation is in agreement 
with the 0(p 4 ) CHPT results in Ref. j2| and can be obtained from a more careful 
discussion jB]. Again this behavior is expected to hold for Im L Tj( 2n ^. For simplicity 
in the following discussion we drop out the spin and isospin indices for partial wave 
amplitudes. 

For the asymptotic behavior of ImRT(s) as s — > 0, in perturbation theory the 
unitarity condition of the partial wave amplitude is satisfied at each order of per- 
turbation expansion on the unitary cut s > 4m\: 

Im R T 2 (s) = , 

n-l 

Im R T 2n (s) = p(s) T 2n -2i(s)T^(s) , for n > 2 . (9) 

i=l 

To be specific, 

Im R T 4 (s) = p(s)T%(s) , 
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Im R T 6 (s) = 2p(s)T 2 (s)Re R T 4 (s) 

= 2p(s)T 2 (s)T i (s)-2tp 2 (s)Tl(s), 



Im R T 8 ( S ) = p(s)(2T 2 (s)Re R T 6 (s) +T 4 (s)T:(s)) 

= 2p(s)(T 2 (s)T 6 (s) + T 4 2 ( S )) - Qtp 2 (s)T 2 (s)T,(s) - Ap\s)T 2 \s) 



n—l n—1 



Im R T 2n (s) = p(s) T 2 n-2i(s)T 2i (s) - 2ip(s) T 2n _ 2i (s)lm R T 2i (s) . (10) 



8=1 1=1 



From above expressions we can easily figure out that as s approaches Im R T 2n (,s) 
is more and more singular when we expand the amplitude to higher orders, because 
of the higher power of p{s). We denote the most divergent term of Im R T 2n (s) 
as IhirT^^s) and the next to leading divergent term as Iim^T^^s). By simple 
deduction, one obtains for n > 3 , 

Im R T«( S ) = (-2 ip (s)T 2 (s)r- 2 p(s)T 2 2 (s), 

lm R T 2 ^(s) = {n-l){-2ip{s)T 2 {s)) n -*T A {s)p{s)T 2 {s). (11) 



From the fact that T 2 (0) is only a nonzero real constant and IhilT 4 (,s) ~ 0((y/— s) 3 ) 
as s —>■ 0_, we conclude that the behaviors of ImLlm R T 2 ^(s) and ImLlm R T 2 ^(s) as 
s — > 0_ are, 

lm L lm R Tt\s) ~ J f^V' ' ■ as .s 0_ if u odd 



Im L Im R T 2 ( J(s) ~ 



R^1m K T$(s) ~ <! / ! \n-i 



Re L Im R T 2 ( J(s) 




as s — ► 0_ if n even, 
0_ if n odd, 

0_ if n even, 
as s — > 0_ if n odd, 
as s — > 0_ if n even, 

0_ if n odd, 

(^7=5 J , as s — > 0_ if n even. 



(12) 

To summarize, the most divergent term of ImLlm R T 2ra (s) and ReLlm R T 2n (s) are : 

, as s — > 0_ if n odd, 

Im L Im R T 2n (s) 

0_ if n even. 
0_ if n odd, 

Re L Im R T 2n (s) ~ { ) v ~ s { n ^ (13) 

0_ if n even. 

From Eqs. (JBJ) and IniL-F and IniL-F contain the same order of divergence 
as p(s)ImLlm R T(s) and ReLlm R T(s), respectively, as s — > 0. For example, the 
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most divergent term of Ihil-F at 0(p 6 ) behaves like 0((l/y/—s) 3 ) and Ihil-F ~ 
0((l/\/— s)), when s — > 0_ and higher order divergences appear as higher order 
perturbation expansions are used. 

It is not difficult to demonstrate that the singular behavior of Ihil-F and Ihil-F at 
s = as discussed above is just a spurious one inherited from the use of perturbation 
theory. 2 To understand this we notice that, for the complete, non-perturbative 
amplitude we have 

which is obtainable from the relation T*(s + ie) = T n (s + ie) = T/S and the 
single channel unitarity relation. As s — > we can obtain from above equation that 
Im R T(s) should not be singular provided that T(s) is not singular, or more precisely, 
Re R T(0) = Tm R T(0) = T(0)/2. Therefore, by simple deduction using Eq. © and 
IullT ~ (^/^s) 3 , one can find that lm L F(s) is 0{{\/ -J^s)) and Im L F(s) is 0{yf^s) 
as s — » 0_. The left hand cut integrals are therefore well defined and are finite except 
at s = for F. 

The Eq. f!14|) not only indicates the correct asymptotic behavior of each quantity 
as s — > 0, but also the reason why the naive use of perturbation results leads the 
left cut integrals divergent - if we could sum up the perturbation series to all orders 
the divergence problem would have disappeared: 



Im R T = J2 Im R T 2n 

n=2 
oo n— 1 



n=2 i=l 

oo oo 

P(*)E E r 2n _ 2i (5)r*( s ) 

i=l n=i+l 

oo oo 

p(«)EE T 2n( S )T 2 *( S ) 
i=l n=l 

/ on \ 2 



= E T *W -2zp(s)EE T 2n( 5 )Im R T 2j ( S ) 

\i=l / »=1 n=l 

(oo \ 2 / oo \ 

^T 2l ( S )J - 2zp( S ) ^E T 2«( s )J Im R T . (15) 

from which we can deduce 

I™ T T^, P( S ) (Ei^l ^2i(s)) 2 , s 

Im R T = ^ Im R T 2n = . 16 

Notice that in deriving Eq. (jlfij) we do not need to use the knowledge on analytic 
continuation. Comparing with Eq. (|14|). the above result is obtained by the simple 

2 Note that this problem will occur not only in ChPT but also in all perturbation theories 
because Eq. JHJ is correct in all perturbation theories and all deductions above are based on this 
equation. In Pade approximations this problem disappears. 
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substitution of T expanded to all orders in perturbation theory into Eq. (JHJ). As 
we discussed earlier, the chiral expansion of T(0) is assumed to be well-defined since 
each T 2n (0) is finite and ~ 0((m^/47r/^) n ). Therefore Eq. (fTflj) suggests the correct 
way of extracting Ihil-F and Ihil-F near s = 0, from a finite order perturbative 
calculation. That is to use 



instead of Y^l=2 I m R^~2i- The former expression removes the spurious divergence in 
the latter introduced by the kinematic factor at s = 0. The Eq. (JT7jl works well in 
the vicinity of s = 0, in other places it should be considered equally well or equally 
bad as the naive expression of perturbation expansion. In our early studies we made 
the naive use of 0{p A ) CHPT to estimate IuilF^') which gives an incorrect Q( t=t ) 



behavior near s' = even though the left hand integral in Eq. is still definable 
except at s = 0. However, the numerical influence to the cut integral at s > 4m^ 
is very small since the integral interval overwhelmed by the spurious divergence is 
very small. Hence numerical results in estimating the left hand integrals are only 
affected very little by the naive use of perturbation theory. 

In Ref. we have discussed the constraint of the Adler zero condition on the 
mi scattering dispersion relation. There we approximately fix the Adler zero at 
s = ml/2 in the I=J=0 channel and force the partial wave S matrix being equal to 
1 at s — rn^/2. However the method is not good enough since, first of all, the Adler 
zero position for the partial wave amplitude is not exactly located at s = rn^/2, and 
secondly, the Adler zero condition actually impose more constraints than what is 
previously considered. 

In partial wave amplitudes, the position of Adler zero can not be exactly given 
because the Adler zero is defined at s = u = t = m\ in the full amplitude and 
after the partial wave projection one cannot fix its exact position. Taking I=J=0 
channel for example, one can find a zero at s — m\j2 in the tree level amplitude 
T% = 2 ™ ^2 ■ If the perturbative amplitude is a good approximation (i.e., converges 
rapidly) to the real situation in the vicinity of s = m\j2 then the Adler zero position 
for the partial wave amplitude, sa, may exist and may be determined by solving 
the equation = T(sa) — T 2 (sa) + T^sa) + ••• using the iteration method: sa — 
— 167r/^(T 4 (m^/2) + ...). For the given perturbative amplitudes from CHPT 
the Adler zero position can be estimated numerically: sa = 0.419 ± 0.058 in unit 
of m\. The error bar appeared in the estimate reflects the uncertainties of coupling 
constants of the chiral Lagrangian [7] . So in the following we will not fix the position 
of Adler zero and instead we use it as a parameter in our fit procedure. 

In Ref. we discussed the role of the Adler zero condition in the global fit. What 
we did is to enforce S in Eq. (9) of Ref. pQ being unity at the zero. Which, however, 
did not make the full use of the Adler zero condition. According to Eqs. (JSJ), (jTljl 
and ReRT(s) = T(l + S)/(2S) it is easy to realize that 



Pis) [Y.UT 2l {s)f 



(17) 



l + 2ip{s) 



F(s A ) = 0, F(s A ) = 1, F'(s A ) = , 



(18) 
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which therefore impose three constraints on the tttt scattering dispersion relations. 
The Eq. f|18|) is helpful in the phenomenological discussions. For example, it can 
be used to make one more subtraction to the dispersion integrals in Eq. and 
two more subtractions to the dispersion integrals in Eq. (J2J). This is appreciable 
since more subtractions can in principle reduce the uncertainties when estimating 
those integrals, which mainly come from high energies. The Eqs. and (J2J) can for 
example be recasted as: 

F(s) = {s-s A )-r— 2 — + (s-4mJ(s-s A )^ 



4m2 - SA K (z 1 / - Ami) (?? ~ * A ) {s - zj 1 ) 

[s-4ml){s- s A ) f ,_ ImF(s') 



7r Jl+r (s' — Aml)(s' — s a )(s' — s) 

Hs) = 1+{s ~ SA)2{s ~ Aml) £ m^m 1 - -a) s (V - 4 <x* - #) 

| {s-Aml)(s-s A f r ^ ImF(s') 

7r Jl+r (s' — Aml)(s' — s A ) 2 (s' — s) 

where a[J denotes the scattering length parameter. Here a(] is no longer a free pa- 
rameter since it is determined by F(s A ) = where F is the one originally defined 
in Eq. (1). 

Using the improved dispersion relations as described above we repeat the fit 
made in Ref. pQ. The fit procedure is the same as before (for example, here we also 
take e = 0.02 which constrains the violation of unitarity ) except that s A is no longer 
held fixed. Taking into account all the uncertainties and variations of parameters 
we arrive at the following results, 3 

M a = 483 ± 13MeV , T a = 705 ± 50MeV ; 

a° = 0.223 ± 0.006 , s A ~ (0.268 - 0.309)m^. (20) 

Note that the scattering length parameter is now in excellent agreement with the 
result of Ref. [HI HH • The above results should replace those given in Ref. P . Notice 
that the error bars given above only represent the uncertainties from our theoretical 
input and does not have a statistical meaning. The position of the Adler zero is 
about 0.3m^ according to our fit, which is not very far from the result from chiral 
perturbation theory. The results on / (980) are also similar to those previously 
obtained. It is worth noticing that the new results are very close to the results 
given in Table 2 of Ref. [2] where the scattering length is constrained by hand using 
the result of Ref. [S]. It was found that the a pole position is rather sensitive to 
the scattering length parameter the correct use of the Adler zero condition is 
therefore crucial in obtaining the correct scattering length parameter and the a pole 
position. 

in a previous version, the fit was performed by inappropriately taking a{] being free, and the 
results were given in Phys. Lett. B549(2002)362 (Erratum). The fit results are nevertheless not 
very sensitive to the different treatment of the ajj parameter. 



7 



We would like to thank Guang-You Qin for helpful discussions. This work is 
supported in part by China National Nature Science Foundation under grant number 
10047003 and 10055003. 



References 

[1] J. Y. He, Z. G. Xiao and H. Q. Zheng, Phys. Lett. B536(2002)59; Erratum 
i.b.i.d. B549(2002)362. 

[2] Z. G. Xiao a nd H. Q. Zheng, Nucl. Phys. A695(2001) 273; Z. G. Xiao and 
H. Q. Zheng, lhep-ph/0103042l 

[3] G. Y. Qin et al, Phys. Lett. B542(2002)89; Q. Ang et al, Commun. Theor. 
Phys. 36 (2000)563. 

[4] J. Y. He, Master thesis, Peking University 2002. 

[5] G. Ecker et al, Phys. Lett. B223(1989)425. 

[6] A. Martin and F. Cheung, Analyticity Properties and Bounds of the Scattering 
Amplitudes, Gordon and Breach, 1969, New York. 

[7] J. Bijnens et al, Nucl. Phys. B508(1997) 263. 

[8] S. Pislak et a/., The BNL-E865 Collaboration, Phys. Rev. Lett. 
87(2001)221801. 

[9] G. Colangelo, J. Gasser and H. Leutwyler, Nucl. Phys. B603(2001)125. 



8 



